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摘 要 : 大 量 观测 数据 分 析 表 明 ,全 球 气 候 正 逐 步 变 暖 。 植 物 物候 现象 是 全 球 自 然 环境 变化 的 指示 
器 。 物 候 对 气候 变化 的 响应 是 研究 全 球 气 候 变 化 的 重要 手段 之 一 。 和 森林 是 全 球 生态 系统 的 重要 组 


成 部 分 ,森林 物候 特征 变化 是 反映 气候 变化 对 森林 生长 影响 的 综合 性 生物 指标 。 利 用 2001 一 2018 
年 MOD09Al 卫星 数据 重建 了 秦岭 地 区 增强 型 植被 指数 (下 朵 7) 序列 ,采用 最 大 变化 速率 和 国 值 法 结 
合 提取 了 秦岭 森林 物候 参数 ,结果 表明 :(1) Whittaker 滤波 法 在 灌木 .农田 森林 3 种 生态 样 地 重 
建 中 稳定 性 较 好 ,在 秦岭 山地 有 较 好 的 适用 性 。(2) 秦岭 地 区 物候 多 年 均值 分 布 同 秦岭 地 区 水 热 
条 件 密 切 ,由 高 海拔 高 山区 到 农耕 区 ,生长 季 始 期 (Start of Growth Season ,SOG) 逐渐 提前 ,生长 季 末 
期 (End of Growth Season, EOG) 2 i 465 , Æ K K (Length of Growth Season, LOG) 由 高 海拔 区 
向 低 海 拔 区 逐渐 加 长 。 秦 岭 浅 山区 和 东部 伏 牛 山 农 耕 带 生长 季 (SOG) 开 始 较 早 ,出 现在 3 ALA, 
高 山区 针 叶 林 带 生长 季 开 始 的 较 晚 ,出 现在 5 月 上 旬 到 中 旬 (120 ~135 d) 之 间 。 生 长 季 末 期 
(EOG) 集 中 出 现在 10 月 ~11 月 初 (270 ~310 d), 高 山区 针 叶 林 带 生长 季 结 束 较 早 , 浅 山 区 植被 生 
长 季 结 束 较 晚 ,普遍 出 现在 11 月 (300 d) 以 后 。 生 长 季 长 度 (LOG ) 分 布 在 150 ~ 200 d 之 间 , 低 海 
拔 地 区 LOG 较 长 ,大 于 180 d, 高 海拔 林 区 生长 季 较 短 LOG 集中 在 150 ~170 d。(3) 年 际 变化 特 
征 :2001 一 2018 年 生长 季 始 期 (SOG) 呈现 提前 趋势 ,其 中 高 海拔 区 提前 明显 ,南北 翡 海拔 低 于 500 
m 部 分 区 域 和 东部 伏 牛 山 少 部 分 区 域 出 现 推迟 。 生 长 季 末 期 (EOG) 呈 现 推迟 趋势 ,其 中 秦岭 北 项 
和 东部 中 低 海 拔 区 域 推迟 明显 ,生长 期 长 度 (LOG) 总 体 呈 延长 趋势 。(4) 秦岭 地 区 近 17 a 气温 呈 


现 上 升 趋势 ,变化 率 为 0.02 ‘CC .a ,降水 呈现 不 明显 的 上 升 趋势 ,日 照 时 数 贝 


| 呈现 明显 的 下 降 趋 


势 , 变 化 率 为 14.6h.a”"。(5) 秦岭 地 区 物候 参数 同 0 C5 CH10 CARIBE AKA RAR 
相关 性 分 析 表 明 , 全 球 变化 下 的 升温 作用 是 影响 秦岭 森林 物候 变化 的 主要 因子 ,升温 作用 导致 
SOG 提前 ,EOG 推迟 、LOG 延长 ,主要 集中 在 秦岭 南北 区 1 000 ~2 000 m 之 间 , 秦 岭 东 部 伏 牛 山 低 


海拔 区 境内 相关 性 最 低 ,表明 受 温度 制约 较 小 。 


X 键 词 : 秦岭 ; 物候 ; 生长 季 末 期 (EOC) ; 生长 季 始 期 (SOG); 生长 季 长 度 (L0G) 


物候 是 指 生物 长 期 适应 气候 条 件 的 周期 性 变 
化 ,形成 与 此 相 适 应 的 生长 发 育 节 律 , 是 反映 生态 系 
统 对 气候 变化 响应 的 最 直观 的 指示 咒 , 通 过 对 和 森林 
物候 的 研究 ,反映 气候 变化 对 树木 生理 、 树 种 组 成 、 
森林 生产 力 等 方面 的 影响 ,对 于 生物 多 样 性 保护 具 
有 重要 的 意义 。 物 候 研 究 是 地 理学 研究 的 热点 之 
一 ,大 尺度 的 植被 生长 季 监 测 是 生态 环境 研究 中 一 
个 重要 科学 问题 "“”。 常 用 的 物候 监测 方法 有 地 面 
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观测 法 .遥感 监测 法 ` 模 型 监测 法 等 。 地 面 观测 法 对 
于 小 区 域 范围 的 物候 研究 具有 时 间 精 度 高 `. 易 于 操 
作 等 优点 ,但 是 难以 反映 区 域 太 度 的 物候 状况 ,卫星 
可 以 在 像 元 尺度 上 监测 植被 群落 的 时 序 特征 , 结合 
地 面 数 据 对 其 参数 化 ,可 以 反映 空间 尺度 物候 特 
1E? A 1990 年 LLOYD 等 "提出 利用 阐 值 法 
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用 ,物候 已 经 成 为 研究 陆地 生态 系统 有 效 的 桥梁 。 
遥感 获取 物候 参数 方法 上 已 有 大 量 研究 ,例如 ,植被 
类 型 生长 的 NDVI 阔 值 法 、 最 大 变化 斜率 方法 等 。 
同时 遥感 监测 的 物候 信息 与 地 面 台 站 的 实际 观测 数 
据 的 相互 关系 越 来 越 受到 关注 。FISHER 26.) 结合 
地 面 物候 观测 数据 .Landsat 和 MODIS 的 NDVI 数据 
集 指出 不 同 源 的 遥感 监测 结果 可 以 描述 同一 区 域 的 
地 面 物候 。KROSS 等 所 研究 表明 时 间 分 辩 率 为 
10 d 和 15 d 的 NDVI 数据 集 适 合用 于 提取 地 表 植 被 
的 物候 信息 。 大 量 研究 证 实 ,在 过 去 的 几 十 年 , 北 半 
球 植物 生长 季 呈 现 逐 渐 延 长 趋势 ,高 纬度 地 区 春季 
物候 的 提前 导致 森林 碳 汇 增加 m1 ,在 我 国 大 部 分 地 
区 ,地 面 和 卫星 观测 数据 印证 了 物候 同样 的 变化 趋 
ga [15-25] A 

秦岭 横 囊 于 我 国 中 部 ,地 处 我 国 暖 温带 和 亚 热 
带 的 生态 过 渡 带 ,是 中 国 南北 地 理 分 界线 ,有 广义 和 
狭义 之 分 :广义 秦岭 西 接 昆 仑 山 ,起 自 岷 山 以 北 , 向 
东 延 伸 经 由 甘肃 境内 、 陕 西南 部 到 达 河 南 境内 伏 牛 
山 一 带 , 狭 义 的 秦岭 只 限于 陕西 境内 渭河 与 汉江 之 
间 的 山区 。 近 年 来 随 着 城市 化 进程 推进 ,人 类 活动 
加 剧 导 致 部 分 区 域 植被 出 现 退 化 -2 。 而 处 于 海 
te 2 800 m 以 上 的 高 山区 植被 ,因为 人 迹 罕 至 ,受气 
候 变 化 影响 较 大 。 目 前 ,涉及 到 秦岭 不 同 海拔 区 域 
的 植被 物候 响应 已 经 有 部 分 研究 成 果 -” ,科研 工 
作者 在 秦岭 山区 植被 退化 及 其 驱动 因素 .以 及 秦岭 
山脉 物候 的 时 空 变 化 方面 进行 了 分 析 , 虽然 涉及 了 
秦岭 物候 的 时 空 变化 ,但 由 于 时 间 序 列 较 短 或 者 过 
于 宏观 ,对 物候 变化 的 原因 及 对 气候 要 素 的 响应 缺 
乏 相 关 的 分 析 。 本 研究 利用 2001—2018 年 卫星 数 


据 ,提取 秦岭 地 区 不 同 海拔 梯度 物候 变化 ,探讨 了 人 
类 活动 和 气候 变化 双重 作用 下 ,不 同 高 程 下 植被 物候 
的 响应 特征 ,为 秦岭 生态 保护 提供 重要 的 理论 依据 。 


1 研究 区 域 概况 


研究 区 域 为 广义 秦岭 ,包括 甘肃 .陕西 .河南 3 
省 的 43 个 县 中 ,境内 自 西 向 东 分 布 着 小 陀 山 ,太白 
山 、 熊 耳 山 、 华 山 以 及 伏 牛 山 等 ,长 约 700 km ,总 体 
上 同 1 月 份 0 CC 等温 线 走向 一 致 ,年 平均 气温 12 ~ 
14 C ,平均 降水 量 为 700 ~900 mm, Al 1 给 出 了 研 
究 区 域 土地 利用 分 布 情 况 , 林 草 分 布 在 较 高 海拔 地 
区 ,耕地 分 布 在 海拔 较 低 的 山谷 地 区 ,中 部 山地 海拔 
垂直 落差 大 ,形成 不 同 植被 类 型 带 谱 , 自 下 而 上 分 
为 :农耕 带 ( 海 拔 200 ~700 m) ;落叶 冰 叶 林带 (海拔 
700 ~1 300 m) ;落叶 栎 林带 (海拔 1 300 ~ 2 300 
m) ;桦木 林带 (海拔 2 300 ~2 800 m) ; 针 叶 林 带 ( 海 
拔 2 800 ~3 200 m) 以 巴山 冷杉 林 和 太白 红 杉 林 为 
代表 , 针 叶 林 分 布 上 界 代 表 森 林 上 限 ; 高 山 灌 从 草包 
带 ( 海 拔 3 200 m 以 上 ) ,秦岭 南北 由 于 气候 条 件 差 
异 , 分 布 稍 有 差别 1。 为 了 提高 验证 精度 ,结合 高 
分 土地 利用 数据 和 地 面 调查 选 出 9 个 典型 生态 样 
地 ,涵盖 森林 灌木 .农田 等 植被 类 型 。 


2 数据 与 方法 


2.1 植被 指数 计算 

采用 NASA LP DAAC (Land Processes Distribu- 
ted Active Center) 提供 的 2001—2018 年 C5 版 本 
MOD09Al 产 品 ,包含 1 ~7 通 道 地 表 反 射 率 、DOY 
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图 1 研究 区 域 土地 利用 图 
Fig.1 Land use changes of Qinling Mountains 
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(Day Of Year) 数 据 , 时 间 分 辩 率 为 8 d, 空 间 分 辨 率 
为 500 m, 秦岭 位 于 H26V5 和 H27V5 区 ,2001 一 
2018 年 共有 822 幅 影 像 ,统一 采用 SIN GRID 投影 。 
为 了 剔除 非 林 区 影响 ,使 用 了 2012 年 MOD12Q1 + 
地 覆盖 数据 ,研究 区 域内 土地 覆盖 类 型 基本 为 林地 
和 耕地 , 占 总 面积 95% 以 上 。 

利用 16 d 间隔 的 MODIS MOD13Q1 产品 统计 
发 现 , 秦 岭 高 植被 覆盖 区 NDVI 会 对 叶绿素 a 响应 
出 现 饱和 现象 1 ,而 EVI 能 避免 饱和 问题 ,减少 大 
气 和 土壤 背景 影响 ,同时 8 d 比 16 d 时 间 精 度 更 高 。 
利用 IDL 语言 结合 MRT 工具 ,编写 批 处 理 程序 ,对 
H26V5 和 H27V5 区 的 反射 率 和 DOY 波段 进行 拼接 
和 投影 ,采用 公式 1 计算 EVI: 


Pir — Prep 

Punir + © X Prep ~ 1. Sparve +1 
式 中 :DNR XP RED ar 分 别 为 MODIS 传感器 大 气 校正 
后 近 红 外 、 红 光 、 蓝 光波 段 反射 值 。 为 了 进一步 提高 
时 间 精 度 , 引 入 像 元 的 DOY 值 , 在 时 间 序 列 重 构 时 
插值 到 日 值 ,避免 波峰 波 谷 点 在 8 d 前 后 波动 。 
2.2 分 析 方 法 

利用 线性 倾向 估计 进行 变化 趋势 分 析 , 变 量 可 
以 看 作 是 时 间 的 一 元 线性 回归 ,线性 倾向 值 用 最 小 
二 乘法 佑 计 ( 公 式 2)。 


EVI =2. 5 x 


(1) 


(2) 


式 中 :B 为 线性 倾向 值 ;X 为 变量 ;i 为 年 份 ;n 为 计 
算 的 年 次 。 当 B>0 时 ,随时 间 ;的 增加 ,X 呈 上 升 
趋势 , 当 B<0 时 ,随时 间 i 的 增加 ,X 呈 下 降 趋势 。 
B 值 大 小 反映 了 变量 上 升 或 下 降 的 速率 , 即 表 示 上 
升 或 下 降 的 倾向 程度 ,同时 采用 相关 系数 的 统计 检 
验方 法 进行 趋势 显著 性 检验 。 

采用 变异 系数 (C,) 评 价 植 被 时 序 稳定 性 :反映 
单一 像 元 序列 离散 程度 ,Cy 值 越 小 ,表明 波动 越 小 
(公式 3)。 


(3) 


用 相关 系数 (RR) 反 映 重建 后 的 时 序曲 线 保持 原 
始 特征 曲线 特征 的 能 力 ,其 值 越 大 ,保持 原始 曲线 特 
征 的 能 力 越 强 (公式 4)。 


> Gea 


R = n n 
Je -8) 之 O: - 7) 


用 回归 估计 标准 差 (RMSE ) 描述 时 间 序 列 之 间 
的 差异 程度 , 值 越 小 , 拟 合 值 的 代表 性 越 强 (公式 


5)。 
> (EVI, - EVI)? 
RMSE = |= 


n 


(4) 


(5) 


2.3 EVI 时 间 序 列 重 构 

为 了 减少 随机 因素 对 序列 的 影响 ,需要 进一步 
对 序列 重 构 ”1 ,常见 算法 有 :(1) 滤波 平滑 法 ,如 
Savitzky-Golay JË yg 、Whittaker JE YE"! 等 。(2) 
函数 拟 合 法 ,如 非 对 称 高 斯 函数 法 ( Asymmetric 
Gaussian ) 、 双 Logistic 函数 法 ( Double Logistic Func- 
tion-fitting) o (3) 频 域 变换 法 ,如 时 间 序 列 谐 波 分 
析 法 (HANTS) 小波 分 析 。 各 种 方法 对 原始 序列 要 
求 苛刻 程度 不 同 ,具有 不 同 适 用 性 ,尤其 是 秦岭 山 
地 ,数据 噪声 点 多 ,这 里 选用 Whittaker 、 双 Logistic 
(DL) HANTS AG 4 种 方法 重建 序列 ,比较 样 地 保 
真性 和 抗 干 扰 性 择优 使 用 。 

图 2 给 出 4 种 算法 2001—2018 年 灌木 .农田 、 
森林 样 地 相关 性 和 RMSE 值 ,可 以 得 到 , Whittaker 
相关 性 最 优 ,DL 和 AG 次 之 日 相当 ,HANTS 略 差 。 

从 RMSE 值 来 看 ,Whittaker 在 灌木 农田、 森林 
3 种 不 同 植被 覆盖 的 重建 中 与 原始 EVI 之 间 的 差异 
性 最 低 ,表明 Whittaker 拟 合 值 的 代表 性 最 强 ,DL 和 
AG 法 次 之 且 相 当 ,HANTS 略 差 。 

对 比 2014—2016 年 4 种 算法 EVI 时 序曲 线 ,得 
到 :(1) 在 灌木 样 地 ,4 种 方法 拟 合 程度 类 似 ,但 在 
峰值 区 ,Whittaker 明显 比 AG DL 和 Hants 抗 干扰 能 
力 强 , 对 于 噪音 多 的 点 表现 出 更 好 的 拟 合 效果 。 
(2) 在 农田 样 地 :生长 季 初 期 , Hants AG 和 DL 拟 合 
结果 均 偏 高 ,生长 季 中 期 至 生长 高 峰 期 4 种 滤波 方 
法 拟 合 结果 一 致 ,并 与 原始 RV/ 时序 曲线 趋势 相符 。 
除 Whittaker 外 的 3 种 滤波 方法 在 2016 年 峰值 期 存 
在 过 度 拟 合 , 且 呈 单 峰 形态 ,但 Whittaker 拟 合 结 
星 双 峰 形态 ,很 好 保留 了 真 值 。(3 ) 在 森林 样 地 :在 
生长 季 初 期 和 枯黄 期 ,4 种 滤波 方法 的 拟 合 结果 基 
本 一 致 ,但 2014 一 2015 年 在 峰值 期 Hants 滤波 受 品 
音 干 扰 拟 合 值 偏 高 ,2016 年 AG 和 DL 滤波 方法 拟 
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Fig.2 Standard deviation and the correlation of the Whittaker, HANTS and Asymmetric Gaussian method over Qinling Moutains 


合 结果 偏 高 , 而 2014 一 2016 年 整个 生长 期 内 Whit- 
taker 与 原始 EVI 时 序曲 线 趋势 最 为 接近 ,很 好 的 保 
留 了 真 值 。 因 此 ,Whittaker 方法 在 4 种 方法 中 表现 
相对 稳定 ,在 秦岭 山地 有 较 好 的 适用 性 , 同 ATKIN- 
SON 等 ' 司 的 结果 一 致 。 

2.4 物候 参数 提取 和 验证 

常见 的 物候 参数 提取 方法 有 羡 值 法 “1 、 比 值 
法 .最 大 变化 速率 法 "| Logistics 函数 法 “等 ,由 
于 秦岭 海拔 落差 大 , 单 用 阔 值 法 对 于 地 形 效 应 体现 
不 明显 ,这 里 选择 最 大 变化 速率 和 疹 值 法 结合 提取 ， 
首先 用 Whittaker 滤波 重建 2001 一 2018 年 EVI 序列 
日 均值 ,计算 曲线 相 邻 点 的 差生 成 EVI 差 值 序列 , 取 
曲线 上 升 期 和 下 降 期 序列 变化 最 大 值 为 SOG (Start 
of Growth Season ) 和 上 OG( End of Growth Season ) 的 
EVI 阔 值 ,最 后 采用 逐年 日 值 序列 对 照 闵 值 指 标 得 
到 SOG 和 EO0G,EOG 和 SOG 差 得 到 LOG ( Length of 
Growth Season ) 。 

从 中 国 物 候 观测 网 获取 了 1963 一 2008 西安 站 
典型 木 本 植物 物候 数据 ,SOG 采用 展 叶 盛 期 木 本 植 
物 的 日 均值 ,EOG 对 应 变色 期 , 展 叶 盛 期 普遍 集中 
在 4 月 上 旬 (95 ~ 100 d) ,全 部 变色 期 集中 在 11 月 
上 人 旬 (300 ~310 d) , 同 卫 星 提取 的 基本 一 致 , 变 化 趋 
WBS EEE ATR 等 人 木 本 植物 展 叶 盛 期 
提前 1.5 d+ a” ,变色 期 延迟 1 d+ a , 同 卫星 提取 
的 变化 趋势 基本 一 致 。 


3 结果 与 分 析 


3.1 2001—2018 年 秦岭 物候 多 年 均值 空间 分 布 
特征 
图 3a ~3c 分 别 给 出 了 2001—2018 年 秦岭 地 区 


RE (LOG ) 多 年 均值 分 布 特征 。 由 图 3, 沿 着 秦岭 主 
要 山脊 由 高 海拔 的 高 山区 到 农耕 区 ,SOG 逐渐 提 
前 ,EOG 逐渐 推迟 ,LOG 由 高 海拔 区 向 低 海拔 区 和 逐 
渐 加 长 ,与 秦岭 地 区 地 形 分 布 较为 一 致 。 秦 岭南 北 
浅 山区 和 东部 伏 牛 山 农耕 带 生长 季 开始 较 早 ,出 现 
在 3 月 上 名 ,高 山区 针 叶 林 带 生长 季 开始 的 较 晚 ,出 
现在 5 月 上 旬 到 中 旬 (120 ~135 d) 之 间 。 生 长 季 末 
期 (E0G ) 集 中 出 现在 10 月 到 11 月 初 (270 ~ 310 
d) ,高 山区 针 叶 林 带 生 长 季 结 束 较 早 , 秦 岭 处 于 浅 
山区 生长 季 结 束 较 晚 ,普遍 出 现在 11 月 (300 d) 以 
Jn. LOG 分 布 在 150 ~200 d 之 间 , 低 海 拔 地 区 LOG 
BEE ,大 于 180 d ,高 海拔 林 区 生长 季 较 短 LOG 集中 
在 150 ~170 d。 随 着 海拔 升 高 SOG 和 LOG 变化 明 
显 ,秦岭 北 坡 浅 山区 带 及 山区 居民 聚集 地 周边 受 人 
类 活动 影响 较 大 ,物候 随 海拔 变化 规律 不 明显 。 
3.2 ”物候 的 年 际 变化 特征 分 析 

图 4a ~4e 分 别 给 出 2001—2018 年 秦岭 地 区 生 
长 季 始 期 (SO0G) .生长 季 末 期 (EOGC) 和 生长 季 长 度 
(LOG) 年 际 变化 趋势 的 空间 分 布 特征 ,图 5 给 出 了 
秦岭 地 区 生长 季 始 期 (SOG ) 趋势 变化 显著 性 检验 分 
布 特征 。 由 图 4:(1) 秦岭 地 区 SOG 分 布 呈 现 提前 
趋势 , 其 中 有 23.45% 象 元 通过 了 显著 性 检验 (P < 
0.01)( 图 5) ,呈现 显著 地 提前 趋势 。 对 于 研究 区 域 
内 的 SOG, 有 95% 的 区 域 提前 1 ~2 d .a1! ,高 海拔 
区 域 SOG 提前 较 明显 ,普遍 在 2 d e a 左右 ,南北 
苞 低 海拔 区 和 东部 伏 牛 山中 部 SOG 呈现 推迟 趋势 ， 
其 推迟 区 域 范 围 同 图 3 中 所 提 及 的 人 类 活动 对 植被 
干扰 较 大 区 域 基 本 一 致 。(2) 秦岭 地 区 EOG 呈现 
推迟 趋势 (图 4b) ,其 中 有 21.45% 的 像 元 通过 了 显 
著 性 为 0.01 的 检验 (P<0.01) (图 略 ) ,呈现 出 显著 
地 推迟 趋势 。 秦 岭 北 草 和 东部 中 低 海 拔 区 域 推迟 明 


生长 季 始 期 (SOG) .生长 季 末 期 (EOG) 和 生长 季 长 


i , 推 后 2 ~3 da (3) 秦 岭 地 区 LOG 变 化 趋 
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图 3 秦岭 森林 物候 多 年 均值 空间 分 布 (2001 一 2018 年 ) (a) SOG;(b) EOG; (c) LOG 


Fig.3 Spatial distribution of forest average phenology over Qinling Moutains from 2001 to 2018 (a) SOG; (b) EOG; (c) LOG 


势 表明 (图 4c) ,2001—2018 年 秦岭 林 区 LOG 总 体 
呈 延 长 的 趋势 ,其 中 70% 的 区 域 延 长 1~3 d.a, 
呈 延 长 趋势 的 区 域 中 有 51% 的 像 元 通过 了 显著 性 
为 0.01 的 显著 性 检验 (P<0.01)( 图 略 ) ,主要 分 布 
于 中 东部 中 高 海拔 区 ,秦岭 地 区 有 30% 的 生长 季 长 
度 缩短 , 其 缩短 日 数 为 -1 ~0 d* a … ,主要 分 布 在 
易 受 人 类 活动 干扰 的 低 海拔 区 。 


4 秦岭 森林 物候 与 水 热 条 件 相 关 分 析 


41 秦岭 地 区 气候 背景 特征 

关于 秦岭 地 区 的 水 热 条 件 变 化 已 经 有 部 分 研 
究 , 主要 集中 在 秦岭 山南 北 苑 春季 和 增 温 效 
应 -3 ， 降 水 呈现 减少 趋势 所 -2 。2000 一 2017 年 
18 a 间 秦岭 地 区 的 的 平均 气温 旦 上 升 趋势 ,变化 率 
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图 4 秦岭 森林 物 变 化 趋势 空间 分 布 (2001—2018 年 ) (a) :SOG;(b) :EOG;(c) :LOG 
Fig.4 Spatial distribution of the trend of forest phenology over Qinling Moutains from 2001 to 2018 (a) :SOG; (b):EOG; (c) :LOC 
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图 5 秦岭 森林 SOG 趋势 变化 显著 性 检验 分 布 (2001 一 2018 年 ) 
Fig.5 Distribution of significance testing of the trend of SOG over Qinling Mountains from 2001 to 2018 


tó wy 


为 0.02 C -a ,平均 气温 最 高 值 为 14.3 C (2013 
年 ) ,最 低 值 为 13.0 °C (2003 年 和 2011 年 )( 图 6) 。 
降水 量 呈 现 出 微弱 的 上 升 趋势 ,变化 率 为 1.5 mm 
a ,降水 量 最 高 值 为 918.0 mm (2003 年 ) ,最 低 值 
为 539.5 mm(2001 年 ) (图 7)。 日 照 时 数 呈现 出 下 
降 趋 势 ,变化 率 为 14.6 h . a” ,日 照 时 数 最 高 值 为 
2 048.4 h (2004 年 ) ,最 低 值 为 1 378.3 h (2016 
年 ) (图 8)。 


14.5 「 年 平均 温度 y=0.021 2x+13.383 
R2=0.0898 
14 | 
p 
z 13.5 上 
> 13 
R 
Br 12.5 
12 
=] N + No} ioe) © N + oO 
[=] © oS © 一 = = = 
© © © © © © © (= =] 
N N N N N N N N N 
年 份 


图 6 2000—2017 年 秦岭 地 区 年 平均 温度 变化 特征 
Fig.6 Annual average temperature over Qinling 


Mountains during 2000 一 2017 
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图 7 2000—2017 年 秦岭 地 区 年 总 降水 量变 化 特征 
Fig.7 Total annual precipitation over Qinling 


Mountains during 2000—2017 
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图 8 2000—2017 年 秦岭 地 区 年 总 日 照 时 数 变化 特征 


Fig.8 Total annual sunshine duration over Qinling 


Mountains during 2000—2017 


4.2 ”秦岭 气候 特征 与 物候 的 关系 
植物 物候 改变 主要 受 环 境 和 气候 变化 影响 , 界 
限 温度 是 指 某 些 重要 物候 现象 或 农事 活动 之 开始 、 
终止 或 转折 点 的 日 平均 温度 ,该 温度 对 植物 的 生长 
有 着 重要 的 意义 。 与 植物 物候 相关 的 几 个 重要 界限 
温度 分 别 为 :0 C ,稳定 大 于 0 °C 的 时 期 为 适宜 农耕 
期 ,其 初 日 与 终日 和 土壤 结 冻 与 解冻 相近 ;5 %C , 稳 
定 大 于 5 的 时 期 为 越冬 作物 生长 活动 期 和 喜 凉 早 
春 作物 的 播种 期 ;10 °C ,稳定 大 于 10 °C 的 时 期 为 越 
冬 作 物 生长 活跃 期 和 喜 温 作物 生长 活动 期 。 为 了 进 
一 步 明确 秦岭 地 区 的 物候 变化 特征 与 界限 温度 的 相 
关 关 系 , 本 研究 分 别 讨论 0% 和 10 所 界限 温度 与 
2001—2017 年 秦岭 物候 变化 特征 之 间 的 关系 。 

计算 了 秦岭 地 区 55 个 站 点 2001 一 2017 年 稳定 
通过 0 C5 CHILO 5 的 积温 .降水 和 日 照 时 数 , 同 
时 提取 了 站 点 周围 5 km 范围 内 的 SOG, EOG, LOG 
均值 ,将 SOG, EOG, LOG 分 别 与 积温 、 降 水 和 日 照 
时 数 做 相关 性 分 析 。 

结果 表明 ,稳定 通过 0 %C 的 积温 .降水 和 日 照 时 
数 与 SOG EOG LOG 的 相关 性 并 不 明显 ,而 稳定 通 
过 5 CA 10 5 的 积温 、 降 水 与 SOG 成 负 相 关 关 系 ， 
与 EOG LOG 呈正 相关 关系 ,稳定 通过 5 CC 和 10 C 
的 日 照 时 数 与 SOG 呈正 相关 ,与 EOG, LOG 成 负 相 
关 关 系 。SOG EOG LOG 与 积温 的 相关 性 要 强 于 降 
水 和 日 照 ,3 个 气候 因子 对 物候 参数 的 响应 顺序 为 : 
积温 > 降水 > 日 照 。 表 明 全 球 变化 下 的 升温 作用 是 
影响 秦岭 森林 物候 变化 的 主要 因子 。 

SOG 总 体 同 稳定 通过 的 10 °C 的 积温 呈现 负 相 
关 关 系 , 其 中 有 13 个 站 点 通过 0. 05 显著 性 检验 
(P <0.05) ,主要 分 布 在 秦岭 山北 芳 的 千 阳 、 凤 翔 、 
收 山 扶风 、 华 阴 、 秦 岭 山 南 世 的 勉 县 、 城 固 、 洋 县 \ 石 
果 、 镇 安 、 商 州 . 丹 凤 、 商 南 等 海拔 高 度 在 1 000 ~ 
2 000 m 之 间 年 均 积温 在 4 000 ~4 500 CHK BK, E 
该 区 域内 ,升温 作用 导致 SOG 提前 。 而 对 于 年 均 积 
温 小 于 4 000 的 高 山区 ,以 及 年 均 积温 大 于 4 800 
C 人 类 活动 影响 较 大 的 区 域 ,SOG 与 积温 相关 性 较 
小 ,表明 在 该 区 域内 , 增 温 对 SOG 影响 不 明显 。 

EOG 和 LOG 总 体 同 稳定 通过 10 °C 的 积温 呈正 
相关 ,其 中 LOG 同 积温 相关 性 较 高 ,15 个 站 点 通过 
0.05 显著 性 检验 (P <0.05) ,通过 检验 的 站 点 分 布 
区 域 同 SOG 与 积温 显著 相关 的 区 域 较为 一 致 。 表 
明 秦 岭南 北条 海 拔 1 000 ~2 000 m 之 间 ( 年 均 积温 
4 000 ~4 500 °C ) 物候 对 积温 响应 较 明 显 ;高 山区 


(年 均 积温 小 于 4 000 °C ) 物候 响应 比较 复杂 ,可 能 
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(SOG) ,生长 季 末 期 (EOG) 和 生长 及 长 度 (L0G) 的 


由 于 针 叶 林 带 的 太白 红 杉 、 巴 西 冷杉 对 温度 变化 异 
常 敏 感 ,而 高 山区 活动 积温 增加 缓慢 ,物候 对 活动 积 
温 响应 较 差 ;秦岭 东部 伏 牛 山 境内 站 点 (年 均 积 温 
大 于 4 900 °C ) 同 界限 温度 相关 和 较 低 ,表明 秦岭 东部 
区 域 受 温度 制约 较 小 。 


5 结论 


和 森林 是 全 球 生 态 系统 的 重要 组 成 部 分 ,森林 生 
长 季 是 反映 短期 或 长 期 气候 变化 对 森林 生长 阶段 影 
响 的 综合 性 生物 指标 。 本 研究 利用 NASA LPACC 
(Land Processes Distributed Active Center ) 提供 的 
2001—2018 年 C5 MOD09Al 产品 ,重建 了 秦岭 地 区 
增强 型 植被 指数 (&W) 时 间 序 列 ,并 评估 了 Whittak- 
er DL HANTS AG 4 种 算法 在 重建 中 的 优 缺 点 ,使 
用 最 大 变化 速率 和 阔 值 法 结合 提取 了 秦岭 物候 参 
数 ,分 析 了 秦岭 地 区 物候 参数 时 空 分 布 和 变化 趋势 ， 
并 讨论 了 秦岭 地 区 物候 变化 与 秦岭 地 区 水 热 条 件 的 
RAR. BRA 

(1) 通过 对 增强 型 植被 指数 EVI 重建 算法 比较 
得 到 :对 于 秦岭 山脉 地 区 ,Whittaker 滤波 法 在 灌木 、 
农田 、 森 林 3 种 不 同样 地 重建 中 相关 性 最 优 ,DL 和 
AG 次 之 且 相 当 ,HANTS MÆ. Whittaker 滤波 法 在 
灌木 农田、 和 森林 3 种 不 同 植被 覆盖 的 重建 中 与 原始 
EVI 之 间 的 差异 性 最 低 , 拟 合 值 代表 性 强 于 DL, AG 
和 HANTS 法 。 此 ,Whittaker 滤波 在 秦岭 山地 有 
较 好 的 适用 性 。 

(2) 2001 一 2018 年 秦岭 地 区 物候 多 年 均值 分 
布 特征 表明 , 沿 着 秦岭 主要 山 痊 由 高 海拔 的 高 山区 


年 际 变化 趋势 表明 ,生长 季 始 期 (SOG ) 分 布 总 体 呈 
现 提 前 趋势 ,有 95% 的 区 域 提前 1 ~2 d- a ' ,高 海 
拔 区 域 SOG 提前 较 明 显 , 普 遍 在 2 d e a EA, A 
北 苞 海拔 低 于 500 m 部 分 区 域 和 东部 伏 牛 山 少 部 分 
区 域 则 呈现 推迟 趋势 ;生长 季 末 期 (EOG ) 总体 呈现 
推迟 趋势 ,其 中 秦岭 北 划 和 东部 中 低 海 拔 区 域 推迟 
明显 , 推 后 2 ~3 d+ a ,生长 期 (LOG) 总 体 呈 延长 
的 趋势 ,其 中 72% 的 区 域 延 长 1~3 da 。 

(4) 秦岭 地 区 55 个 站 点 2000 年 到 2017 年 平 
均 温 度 呈 现 上 升 趋势 , 增 温 幅 度 为 0.02 °C + a! ,年 
降水 总 量 呈 现 不 明显 的 增加 ,而 年 总 日 照 时 数 则 呈 
现 明显 的 下 降 趋势 , 下 降幅 度 为 14.6h a`, 
(5) 水 热 条 件 的 变化 ,进一步 影响 了 秦岭 地 区 
物候 的 变化 特征 。 秦 岭 地 区 55 个 站 点 的 稳定 通过 
0 C5 JC 以 及 稳定 通过 10 °C 的 积温 .降水 和 日 照 时 
数 与 秦岭 地 区 物候 参数 相关 关系 分 析 表 明 ,稳定 通 
过 0 5 的 积温 、 降 水 和 日 照 时 数 与 SOG EOG, LOG 
的 相关 性 并 不 明显 ,而 稳定 通过 5 C A 10 C 的 积 
温 、 降 水 与 SOG 成 负 相 关 关 系 ,与 EOG LOG 呈正 相 
关 关 系 , 稳 定 通过 5% 和 10% 的 日 照 时 数 与 SOG Œ 
正 相 关 , 4G EOG, LOG 成 负 相 关 关 系 。SOG 、EOG、 
LOG 与 积温 的 相关 性 要 强 于 降水 和 日 照 ,3 个 气候 
因子 对 物候 参数 的 响应 顺序 为 :积温 > 降水 > 日 照 。 
升温 作用 导致 SOG 提前 、EOG 推迟 、LOG 延长 , 秦 
岭南 北 芒 海拔 1 000 ~2 000 m 之 间 物 候 对 10 °C FR 
温 响 应 较 明 显 ,高 山区 物候 响应 比较 复杂 ,可 能 由 于 
针 叶 林 带 的 太白 红 杉 、 巴 西 冷 杉 对 温度 变化 异常 敏 
感 ,而 高 山区 活动 积温 增加 缓慢 ,物候 对 活动 积温 响 


到 农耕 区 ,生长 季 始 期 (SOGC ) 逐 渐 提 前 ,生长 季 末期 
(EOG) EMES ,生长 季 长 度 (LOG ) 由 高 海拔 区 向 
低 海 拔 区 逐渐 增长 ,与 秦岭 地 区 地 形 分 布 较为 一 致 。 
秦岭 南北 浅 山 区 和 东部 伏 牛 山 农 耕 带 生长 季 开始 较 
时 ,出 现在 3 月 上 多 ,高 山区 针 叶 林 带 生 长 季 开始 的 
较 晚 ,出 现在 5 月 上 旬 到 中 旬 (120 ~ 135 d) 之 间 。 
EOG 集中 出 现在 270 ~310 d, 高 山区 针 叶 林 带 生 长 
季 结 束 较 早 ,秦岭 处 于 浅 山 区 生长 季 结 束 较 晚 , 普 i 
出 现在 11 月 (300 qd) 以 后 。 生 长 季 长 度 (LOG ) 主要 
分 布 在 150 ~ 200 d 之 间 , 低 海拔 地 区 LOG 较 长 ,大 
于 180 d ,高 海拔 林 区 生长 季 较 短 LOG 集中 在 150 ~ 
170 d, 

(3) 2001 一 2018 年 秦岭 地 区 生长 季节 始 期 


应 较 差 ,秦岭 东部 伏 牛 山 境 内 站 点 与 界限 温度 相关 
性 较 低 ,表明 秦岭 东部 区 域 受 温度 制约 较 小 。 
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Abstract: Based on the product of C5 MODO9A1 from Land Processes Distributed Active Center ( NASA 
LPACC) from 2001 to 2018, the Enhanced Vegetation Index (EVI) over Qinling Mountains, the geographical 
boundary of south-north China, is reconstructed. The four reconstruction methods , Whittaker, DL, HANTS and AG 
algorithms were evaluated. Then the phenological parameters of Qinling Mountains was extracted with the maximum 
change rate and threshold method. We got the followed results. The Whittaker filter method is generally the best a- 
mong the four methods in reconstructing the EVI of shrubs, fields and forests over Qinling Mountains. And the EVI 
from the Whittaker filter method shows the lowest difference from the original EVI in the reconstruction of shrub, 
farmland and forest cover with higher representative than that from the DL, AG and HANTS algorithms. Along the 
main ridge of the Qinling Mountains from the high elevations to the farming area, the start of growth season is gradu- 
ally moved earlier ,the end of growth season is gradually delayed and the length of growth seasons become much lon- 
ger. The growing season in the northern and southern shallow-mountain areas and Funiu Mountains starts in early 
March , while it begins in early and mid-May (120 - 135 d of the year) in the coniferous forest belt. While the end 
of growth season in the coniferous forest belt ends earlier than that of shallow-mountain areas which mainly appear 
in early of November, and both of them is mainly in the range of 270 -310 d of the year. The length of the growth 
(LOG) season over Qinling Mountains is mainly about 150 - 200 d,and the LOG at lower elevation ( about 180 d) 
is longer than that of forest regions at higher elevation which is about 150 - 170 d. The distribution of the initial pe- 
riod of SOG over 95% areas of the Qinling Mountains generally shows an advance of 1 -2 d, especially for the high- 
altitude area, which shows an obvious advance trend with about 2 d - a ,while there is a delay trend over the areas 
with the elevation lower than 500 m and the Funiu mountains. There is a delayed trend of the growth seasons 
(EOG) over Qinling Mountains especially over the north of Qinling Mountains and low-altitude areas (about 2 -3 
d + a`! ). Furthermore, the length of growth seasons over Qinling Mountains become longer, and the change range of 
72% of the forest is about 1 -23 d+ a .There is a slightly increased trend of average temperature with a rate of 
0.02 °C + a ',while the total annual sunshine duration is decreased with a rate of 14.6 h + a ' over the Qinling 
Mountains. The change of hydrothermal conditions further affected the phenological parameters of Qinling 
Mountains. The correlation analysis between accumulated temperature , total precipitation and total sunshine duration 
and the phenological parameters at the temperatures of 0 °C ,5 °C and 10 C respectively , which shows that the re- 
sponse sequence of the phenological parameters to the three climatic factors is; accumulated temperature > precipi- 
tation > sunshine. The increasing of temperature slightly leads to an earlier SOG ,a delayed EOG and a longer LOG, 
especially in the areas with the elevation between 1 000 and 2 000 m over the southern and northern part of Qinling 
Mountains. While for higher mountains the correlations becomes more complex. The response of phenological param- 
eters to accumulated temperature over Funiu Mountains in the east of Qinling Mountains is not obvious and there 
may be some other restriction factors. 
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